
m09 = read_mrt('./local/mrt/Marshall2009_Table1.txt')

restore,'./irdc_dist_model/bgps_rb3.sav',/ver

n = n_elements(rb3)

m09ind = lonarr(n)

FOR i=0L, n-1 DO BEGIN
   
   gcirc,2,rb3[i].l,rb3[i].b,m09.glon,m09.glat,dis
   
   d_min = min(dis,mind)
   m09ind[i] = mind
   
   IF d_min LE 60. THEN print,d_min, rb3[i].morph_d,m09[mind].dist
   
ENDFOR




END
